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ABSTRACT 

The observational characteristics of quasi-periodic oscillations (QPOs) from accreting neutron stars strongly 
indicate the oscillatory modes in the innermost regions of accretion disks as a likely source of the QPOs. The 
inner regions of accretion disks around neutron stars can harbor very high frequency modes related to the 
radial epicyclic frequency k. The degeneracy of n with the orbital frequency VL is removed in a non-Keplerian 
boundary or transition zone near the magnetopause between the disk and the compact object. We show, by 
analyzing the global hydrodynamic modes of long wavelength in the boundary layers of viscous accretion 
disks, that the fastest growing mode frequencies are associated with frequency bands around k and k ± 57. The 
maximum growth rates are achieved near the radius where the orbital frequency 51 is maximum. The global 
hydrodynamic parameters such as the surface density profile and the radial drift velocity determine which 
modes of free oscillations will grow at a given particular radius in the boundary layer In accordance with the 
peak separation between kHz QPOs observed in neutron-star sources, the difference frequency between two 
consecutive bands of the fastest growing modes is always related to the spin frequency of the neutron star. This 
is a natural outcome of the boundary condition imposed by the rotating magnetosphere on the boundary region 
of the inner disk. 

Subject headings: accretion, accretion disks — stars: neutron — stars: oscillations — X-rays: stars 



1. INTRODUCTION 

Wave modes in the boundary region in the inner accre- 
tion disk are a likely source of the distinct narrow frequency 
bands of quasi-periodic oscillations (QPOs) (van der Klis 
2000) from neutron star sources in low mass X-ray binaries 
(LMXBs). The dependence of QPO frequencies on accretion 
rate suggests that the observed QPOs are connected with ac- 
cretion disk modes. Psaltis, Belloni & van der Klis (1999) 
showed the existence of correlations between the frequencies 
of different QPO bands extending over a wide span of fre- 
quencies. The same correlation encompasses black hole as 
well as white dwarf and neutron star sources. This strongly 
implies that the frequency bands are determined by the oscil- 
lation modes of the accretion disk. The nature of the compact 
object, whether it is a black hole, white dwarf or neutron star, 
probably plays a role in exciting the same disk modes through 
possibly different mechanisms in different types of compact 
sources. The QPO bands modifying the X-ray luminosity are 
likely to belong to the boundary region between the disk and 
the compact object. For magnetic neutron stars, this bound- 
ary region is shaped by the interaction of the disk with the 
magnetosphere. The region where the rotation deviates from 
Keplerian flow is not necessarily very narrow; this boundary 
zone may have a size as large as a few tenths of the inner 
disk radius. We shall therefore use the terms "boundary" and 
"transition" region interchangeably. 

Initial explorations of disk modes underlying QPO fre- 
quency bands from LMXBs were provided by Alpar et al. 
(1992) and Alpar & Yilmaz (1997). The characteristics and 
possible excitation mechanisms of thin-disk oscillations were 
reviewed and discussed by Kato (2001). In many current 
models of QPOs, especially those of neutron stars, their fre- 
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quencies are identified with test-particle frequencies (Stella, 
Vietri & Morsink 1999; Abramowicz et al. 2003). Most parts 
of an accretion disk, except, significantly, the transition re- 
gions at the inner boundary of the disk, are characterized by 
Keplerian rotation rates to a very good approximation. At 
any radius beyond the transition region, acoustic, magneto- 
acoustic and viscous corrections to the test-particle Keplerian 
orbital frequency are negligible. In the outer disk regime for 
disks around neutron stars and white dwarfs (and far enough 
from a central black hole), test-particle frequencies for radial 
and vertical perturbations of the orbit are degenerate with the 
Keplerian orbital frequency. The degeneracy is removed if 
the effects of general relativity are important. A Newtonian 
field of tidal or higher multi-pole structure would also lead 
to a split in the degeneracy. However, in Newtonian grav- 
ity, distortions of stellar shape even for the most rapidly ro- 
tating neutron stars will not introduce a significant level of 
non-degeneracy between the frequencies of test-particle os- 
cillations that is comparable to the difference between the ob- 
served QPO frequency bands. Models employing the Keple- 
rian frequency have been relatively successful in interpreting 
QPO frequency correlations. The Kepler (test particle) fre- 
quency interpretation has also been employed to place con- 
straints on the masses and spins of compact objects. Empir- 
ically, predictions based on test-particle frequencies applied 
to QPO frequency correlations deviate from observations at a 
level of about 10%. 

The Keplerian frequency represents the basic imprint of ro- 
tation on all the dynamical responses of the disk. Albeit sim- 
ple, the identification of QPOs with test-particle frequencies 
has a number of shortcomings. First, precisely in a bound- 
ary layer where QPOs are expected to be excited, the dynam- 
ical frequencies of oscillations are not the test-particle fre- 
quencies. In the boundary layer, viscous and magnetic forces 
lead to deviation of orbital frequencies from Keplerian test- 
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particle frequencies (see, e.g., Erkut & Alpar 2004). The re- 
sulting band of non-Keplerian rotation frequencies entail vis- 
cous and acoustic response and couplings between adjacent 
rings of the disk fluid. Hydrodynamic effects are therefore es- 
sential for an understanding of the frequency bands character- 
izing the boundary layer. In particular, the degeneracy of test- 
particle frequencies in the non-relativistic regime is lifted in 
the hydrodynamic boundary layer, where the radial epicyclic 
frequency is in fact larger than the orbital frequency. This 
simple observation (Alpar & Psaltis 2005) has important con- 
sequences for the interpretation of kHz QPO frequencies, in 
particular for the constraints on the neutron star mass-radius 
relation derived from kHz QPOs. 

Second, test-particle frequencies do not distinguish be- 
tween azimuthal sidebands. While kinematic models of QPOs 
like the beat-frequency model involve one specific band of 
frequencies, say lo, a large number of sidebands, of frequen- 
cies Lo,„ = ijj- mfl, where ft is the orbital frequency, are al- 
lowed by azimuthal symmetry. Arguments as to why only 
one or two QPO bands are excited must involve choices im- 
posed by the symmetries of the interaction between the mag- 
netosphere and the accretion disk boundary, leading to reso- 
nances with particular frequencies at one or more radial re- 
gions in the disk. Without resonances, test-particle frequency 
spectra present no distinction between the modes. The re- 
alistic hydrodynamic modes may, in some parameter ranges, 
distinguish between both fundamental modes and all their az- 
imuthal sidebands through the different growth or decay rates 
of these modes, even in the case of free oscillations. A reduc- 
tion of the number of relevant (easily excitable) modes of free 
oscillations is certainly an important task for an understand- 
ing of accretion disks around neutron stars, white dwarfs and 
black holes. In the case of black holes, this classification of 
free hydrodynamic modes in terms of their growth rates is 
even more important as resonant excitation by the black hole 
is not available. 

In this paper, we study the global modes of free oscilla- 
tions at some position r in the inner disk-boundary or transi- 
tion region by analyzing the perturbed dynamical equations 
of a hydrodynamic disk, including pressure gradients, vis- 
cous and magnetic stresses. We consider the case of a disk 
around a neutron star with a magnetosphere. The neutron star 
is taken to be a slow rotator, the star's rotation rate is less 
than the value of the Keplerian rotation rate at the inner ra- 
dius of the disk; fJ* < Q,-^{rin). The actual rotation rate of the 
disk is set by the boundary condition f2(rin) = 51*. Thus in the 
steady state solutions for the disk, the rotation rate at the inner 
edge of the disk and in the transition region beyond is sub- 
Keplerian. The innermost disk radius rin determines the disk 
magnetosphere interface which may be subject to magnetic 
Rayleigh-Taylor or interchange instabilities. This instability 
operates mainly at the disk-magnetosphere interface because 
of the sharp density contrast across the radius r = ri,,. The in- 
ner disk region for r > consists of a non-Keplerian bound- 
ary region that joins the outer Keplerian disk with a continu- 
ous density distribution. As we will see in § 13.21 the surface 
density profile throughout the non-Keplerian boundary region 
increases with decreasing radii. Thus, our boundary region is 
not subject to interchange-like instabilities. 

The free oscillation modes we explore by perturbing a 
steady state solution embody information about the stability 
of the disk through their growth or decay rates. Growing 
modes will provide a clue as to the origin of the QPO fre- 
quency bands, which could be excited by resonant forcing of 



the disk by time dependent interactions with the star's mag- 
netosphere. We find, as expected, that growth or decay rates 
are determined by the dynamical effects of viscosity, with 
an additional dependence on the sound speed and the den- 
sity and rotation-rate profiles in the non-Keplerian boundary 
region. For some boundary region parameters, not all funda- 
mental modes of free oscillations and not all azimuthal side- 
bands grow, and among the growing modes the growth rates 
differ For other choices of boundary region parameters the 
growth rates of all sidebands are similar. The frequencies 
of the modes differ from test-particle frequencies by amounts 
that can be several times larger than the corresponding QPO 
width. Most importantly, for a reasonable set of accretion 
disk parameters, we can show that only a certain few of the 
hydrodynamic free oscillation modes will grow, and that the 
frequency bands of these oscillations can be associated with 
the observed frequency bands. The nature of free oscillation 
modes in the boundary region depends on whether the neu- 
tron star is rotating at a rate slower or faster than the rotation 
rates prevailing in the inner boundary of the disk. In this pa- 
per, we will consider the more common and straightforward 
case of slow rotators, the case when the neutron star rotation 
rate il» is less than ilait'm), the Kepler rotation rate at some 
representative inner disk radius ri„. An investigation of forced 
resonant excitation of these prevalent modes, and the compar- 
ison and association with observed QPO bands will follow in 
a subsequent paper. 

§ 2 lays out the basic assumptions and equations, § 3 dis- 
plays the mode analysis for the free oscillations, including 
a discussion of hydrodynamic effects, frequency bands, and 
growth rates. In § 4, we discuss the results and present our 
conclusions. 

2. BASIC EQUATIONS AND ASSUMPTIONS 

We consider the excitation of oscillations in a geometrically 
thin disk in vertical hydrostatic equilibrium. We write the 
continuity and the radial and angular momentum equations 
in cylindrical coordinates, integrated vertically over the disk 
thickness, as 

/9S 1 1 ^ 

— + -— (rSv,)+-— (l]v^) = 0, (1) 
ot r or r ocp 

and 

at or r ocp r Sr ocp 

Here, H is the half-thickness of the disk, v,. and are the 
radial and azimuthal components of the velocity field in the 
disk, $ is the gravitational potential, 

S = / pdz (4) 

J-H 

is the surface mass density, 

n= / Pdz (5) 

J-H 

is the vertically integrated thermal pressure, p is the mass den- 
sity, P is the sum of the gas and radiation pressures, and 

(Fj) ^ ^ fpjpdz (6) 

J —H 
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are the vertically averaged sums of the magnetic and viscous 
forces per unit mass, and the subscript j stands for the r or cj) 
component of the force. 

We only consider the 0-component of the viscous force as 
the dominant shear stress between adjacent layers, because of 
our assumption that the disk is geometrically thin (Shakura & 
Sunyaev 1973). We write the sum of the viscous and large- 
scale magnetic forces (per unit mass) in the radial and az- 
imuthal directions as f,. = f and = F^" + F^'^° , respec- 
tively. 

The steady, axisymmetric, equilibrium state obeys 



27rrI]()Vro = M, 

z=H 



VrQ—^-Q. r 

ar 



dr J So dr 



and 



(7) 

+ {FDo^ (8) 
(9) 



where M is the constant mass-inflow rate, v,o < 0, v^o = ^r, ^l 
is the angular velocity of the disk plasma, and k = [2fl(2r2 + 
rdfl/dr)y^^ is the radial epicyclic frequency. 

We introduce perturbations to this axisymmetric fluid dis- 
tribution and use the subscripts "0" and "1" to denote the equi- 
librium and perturbed states, respectively. The linearized per- 
turbation equations follow from equations ([T])-© as 

^ + - I- [KSoVH + Si VK))]+ - T^(SoV^i) = 0, (10) 
cLt r or r oo 



— 2S 2V0i + -^(VrOVrl ) = — 

dr or So 



Si dU, 







So dr 



dr 



(11) 



and 



dr 2U r Or 



1 dlVi 

Sor d<j) 



(12) 



where d / dr = d / dt+Vld / d(t>\s, the Lagrangian derivative fol- 
lowing the motion of the fluid in the azimuthal direction. 

In the following, we study perturbations varying on 
timescales shorter than the thermal timescale and thus expect 
that there is no energy exchange between adjacent fluid ele- 
ments. In this case, the fluid elements respond adiabatically to 
density and pressure perturbations and, at the same time, the 
disk thickness H, which is determined by vertical hydrostatic 
equilibrium, responds to the change of pressure. In general, 
we can write (Stehle & Spruit 1999) 



n, 



So A d\np) 



(13) 



where F is the index of a polytropic equation of state. The 
terms in the last parenthesis is of order unity. Because we will 
not be calculating explicitly changes in the vertical equilib- 
rium during the oscillations, we will set the last term to unity. 
Therefore, we will be studying effectively isothermal modes. 

We consider long wavelength perturbations to analyze the 
lowest order normal modes extending globally over the radial 
distance scale rout of the disk. For the stability of such global 
modes the effect of cylindrical geometry and of the structure 



of the unperturbed disk state must be taken into account. We 
perform a local mode analysis. Our approach differs from the 
usual local mode analysis of short-wavelength perturbations 
of radial dependence e'*'' in a WKB approximation. Consid- 
eration of modes with radial wavenumbers would be appro- 
priate for the analysis of local disk modes in the short wave- 
length regime. Instead, we study the stability of the global 
disk modes at each particular radius r in the boundary layer In 
this region, the unperturbed disk quantities such as the surface 
density So and the rotation rate Vl change on a length scale 
(Ar)BL < Hn- Since the variation of the global long wave- 
length perturbations has a much longer range roui ^ (Ar)BL 
in the radial direction, we write equations (fT0l)- (fT2l) at any 
particular radius r in dimensionless form, neglecting the r- 
dependence of the perturbations: 



+ ( 1 + /3)aM,- + + a = 0' 

(XT ocp 
dr 



du^ 

dr 2\l 



-a 



da 



(14) 



(15) 



(16) 



Here, a = Si /So is the dimensionless density perturbation, 
M,. = Vri/c, and = v^i/c, are the dimensionless velocity 
perturbations, /^.^ = (iv,^) i / Q are the force perturbations in 
units of frequency, and = -v^o/r and 57, = c jr are the 
typical frequencies associated with the radial accretion and 
the effective sound speed c, = (IIo/So)'''^, respectively. The 
parameter 

represents the surface mass density profile at a particular ra- 
dius r in the inner disk. Note that (3 crucially represents the 
structure of the unperturbed disk. This differs from the usual 
local mode analysis in which the radial variation of the unper- 
turbed background quantities such as So is neglected while 
the radial derivatives of the perturbations are taken into ac- 
count in terms of their radial wavenumbers. 

We look for solutions of the above linear set of equations of 
the form 

qici), t) = Y, f q,niuj)e"-"'^-"'>du;, (18) 

where the integer m>\ represents non-axisymmetric modes 
corresponding to the prograde motion of perturbations in the 
azimuthal direction. The Fourier decomposition of equations 
dl-dlSll yields 

[i{mVl- uj) + I3VL^] a,„ + (l+ l3)flsUr,m + '"^^s'^<t>,m = 0, (19) 

- (3n,a„ + [iimn -uj) + il+P)n,] u„n-2^u^.m = /™s , (20) 



imn,(7„ + ^Mr,m + [i(mn -LU)- Q^] u^^„ = f^'^„, + /™f . (2 1 ) 

To proceed further in the linear mode analysis, we need to 
specify the force perturbations fr,„, and f^j„ in equations ( |20] i 

and (Ell- 
in this paper, we discuss the free oscillation modes. In a 
subsequent paper, we will study oscillations forced by exter- 
nal perturbations, as in the case of a compact object with an 
oblique large-scale magnetic field. 
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3. GLOBAL FREE OSCILLATIONS IN A VISCOUS ACCRETION DISK 

When there are no perturbations introduced by large-scale 
magnetic fields, i.e., when /i™^ = = 0, the only non- 
negligible force perturbation is due to the kinematic viscosity 
V. In the absence of an analytic model for the effective kine- 
matic viscosity in a turbulent shearing flow, we adopt a sim- 
ple damping force prescription, i.e., F^'*^ = -v^/f^, where is 
the timescale for viscous accretion. In order for our prescrip- 
tion to yield the order-of-magnitude estimate for the viscous 
timescale, ^ p- jv ^ rj |vro|, we write 



■■Iv- 



(22) 



where ^i, is a dimensionless factor. Note that, for this pre- 
scription of the viscous force to satisfy equation (|9]l for the 
steady equihbrium state, it is necessary that 

> = — (23) 

To linear order in the perturbed quantities, equation ( |22] | 
yields (^0)) = -^^^^yV^x. The zeroth-order magnetic force 
enters the perturbation equations only through the parameter 
7i/ as given in equation (|23] |. 
For the case of interest, equations (fT9]l-(l2TI) become 

[/(mil - w) + Om + (1 + /3)^sUr,m + imflsU^ ,„ = 0, 

— prisCFm + [iimfl — Uj)+(l + /3)r2 J Mr,m — '^^Ucp.m = 0, 



(24) 
(25) 

(26) 



yielding a dispersion relation of third order in uj. We pro- 
vide the general expressions for the dispersion relation and 
the three complex solutions in the Appendix. 

For illustrative purposes, we write the expressions for the 
three frequencies to leading order in rij/il, i.e., for small hy- 
drodynamic corrections. The solutions of astrophysical inter- 
est may actually be found in situations beyond this approxi- 
mate regime, as we shall discuss in § 3.2. However, the ap- 
proximate solutions shed light on the basic parameters that 
determine the frequencies and growth rates of the modes. In 
the limit of small hydrodynamic corrections, i.e. ei, 62, and 
£3 ^ 1 (see Appendix), the axisymmetric modes (m = 0) have 
the frequencies and growth rates 

=±K±f3(i+i3) \ 2^ la 



Re 



U! 



± 



Im 



2k 



3(7,-/3- 1)-(7,-/3)^ 



6k 



,(0) 



Re 



,(0) 



and 



Im 



,<0) 



= 0, 



n,, (27) 



(28) 



(29) 



(30) 



The corresponding expressions for the non-axisymmetric 
modes, with positive or negative integer values of the az- 
imuthal wavenumber m are 



Re 



(m) 
-'1,2 



± 



3(7,-/3- l)-(7.-/9)' 



^1^^.(31) 



Im 



and 



Im 



(m) 



Re 



,(m) 



(m) 



7. + 2/3 



7. + 2/3 



(32) 



(33) 



(34) 



Modes whose real frequencies differ only by a sign are of 
course identical. We note that in these approximate expres- 
sions to leading order in fl^ the imaginary parts have no de- 
pendence on m. This means that the fundamental modes 
and the infinite sequence of sidebands shifted in frequency 
by ±|m|f7 from the fundamental modes of interest, all have 
the same growth or decay rates in the limit of small hydrody- 
namic corrections, i.e. 51, ^ and fij <C fl- Any distinction 
between the free oscillation modes arises only in a regime in 
which acoustic and viscous hydrodynamic effects are impor- 
tant. 

It is possible to check the consistency of our stability anal- 
ysis using equations (l28b . (l30b . ( |32] |. and ( l34b in the limit of 
non-magnetic local modes for which the radial gradient in the 
background surface density at any given location is negligi- 
ble, that is, dln'Eo/dlnr ^ 1. Note that none of the modes 
grows if the local gradient of the surface density is neglected 



(/3 = 0) since 7, >0 for (F" 



: 0. There is no instability for 



local hydrodynamic modes in the long wavelength limit. The 
global modes grow in the inner disk where the radial gradient 
of Eo cannot be neglected {f3y^0). 

3.1. Hydrodynamic Corrections to Test-particle Frequencies 

In the ideal case of fi, = fij = 0, the three characteristic 
frequencies coincide with appropriate combinations of test- 
particle frequencies in the disk, i.e., wj'"' ~ k+ \m\Q, Wj"' — 
K- \m\n, and Wj"' ~ \m\fl. This justifies, to zeroth order, the 
identification of the observed QPO frequencies with frequen- 
cies of test particles in kinematic models of the QPOs. For 
example, in the sonic-point model (Miller et al. 1998), the 
higher kHz QPO would be the mode at the sonic radius. 
In the relativistic precession model (Stella et al. 1999), the 
two high-frequency QPOs would be the Wj'^ and Wj'^ modes. 

In the absence of external forcing, high-frequency modes 
in an accretion disk can be excited only in the presence of 
viscosity. This is a well-known result, discussed in detail by 
Kato (2001). Because the disk cannot respond thermally at 
the high-frequencies of interest here, viscosity can excite the 
modes only through its dynamical effect (mechanism II in 
Kato 2001) and not through its thermal effect (mechanism I 
in Kato 2001). Thus free oscillation modes in the accretion 
disk modes can have positive growth rates only if il, is not 
zero, i.e., through the effect of viscosity, and, hence, only if 
the hydrodynamic contributions to their frequencies are op- 
erative. Indeed, all models of hydrodynamic accretion disk 
modes show that non-negligible hydrodynamic corrections 
modify the test-particle frequencies of growing modes (see, 
e.g.. Wagoner 1999; Kato 2001; Psaltis & Norman 2000). 

We estimate the expected hydrodynamic corrections to the 
oscillatory mode frequencies typically identified with high- 
frequency QPOs using the scaling of the viscous and thermal 
frequencies in an alpha disk, i.e.. 



(35) 
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and 



H 



n. 



(36) 



The 
Re 



fractional hydrodynamic correction to, e.g., the 
Wj'^ ~ «; - fi mode is determined by the acoustic re- 
sponse provided that a is not close to unity: 

2 /rr\2 _ . /r^ \ 2 



6ul_ 
, ,(1) 



H 



(37) 



The growth rate of the same modes is ~ Q,^ and hence the 
quality factor Q of the corresponding QPO will be at most 



Q' 



-n 



H 



n 



1-^ 
n 



(38) 



and the resulting width of the QPO will be ~ — 
Slj/. Combining equations ( l37b and (|38] ) we obtain for the 
hydrodynamic correction to the mode frequency in units of 
the width of the corresponding QPO the relation 



Suj'. 



(1) 



1 fn 



UJ. 



(1) 



/Q 



(39) 



It is clear from this expression that, even for relatively large 
values of the parameter a ~ 0.1, the QPOs will have frequen- 
cies that are distinct from the corresponding test-particle fre- 
quencies by amounts that are several times larger than the 
QPO widths. 




r /r ■ 

Fig. 1. — The radial profiles of the orbital frequency Q(r), epicyclic fre- 
quency K(r), Keplerian frequency S1k(''). ™d the vertically integrated dy- 
namical viscosity f(r) throughout the inner disk. This example is based on a 
typical boundary region model from Erkut & Alpar (2004). The frequencies 
ai'e given in units of QKCnn)- In this particular example, the actual orbital 
frequency at the inner disk boundary and the rotation rate of the neutron star 
are half the Keplerian value at fi„; Sl(rjn) = = f2K(''in)/2. 



3.2. Excitation of Global Hydrodynamic Free Oscillation 

Modes 

The growth rates of free oscillation modes are determined 
by a number of parameters characterizing the boundary re- 
gion in the inner disk. In this region, the angular velocity of 
the disk matter makes a transition from the Keplerian rotation 
to match the stellar rotation rate. Our analysis depends on 
the local values of these parameters at any particular radius r 



within this transition or boundary region. The main param- 
eters are the ratio of the radial epicyclic frequency k to the 
orbital angular frequency ft, the radial surface-density profile 
/3 (equationfTTTi. il^/ils, the ratio of the radial drift velocity to 
the sound speed, and il^/^, the inverse timescale or the typi- 
cal frequency associated with the sound speed in units of the 
angular velocity. In this section, we relate these parameters to 
the angular velocity profile i7(r) and the dynamical viscosity 
i/Eo in the accretion disk. We then study the local excitation 
of the modes in the inner disk or boundary region where the 
hydrodynamic effects are important. 

The ratio of the epicyclic and orbital frequencies is re- 
lated to the orbital frequency profile i7(r) through the relation 

lt/lnf7\ 

= 4(1 + -^^ I. (40) 



2 d\nr 



In a non-Keplerian boundary-transition region of the accre- 
tion disk around a neutron star that is a slow rotator, Vl{r) 
is less than the Keplerian value 51k('")- Proceeding from the 
Keplerian outer disk through the transition region, Vl{r) goes 
through a maximum and then decreases to match the star's 
rotation rate J7» at the inner radius of the disk. The epicyclic 
frequency k, degenerate with JIk in the outer disk, increases 
through the transition region. The ratio n/Vl equals 2 at the 
radius where VL is maximum, has values between 1 and 2 in 
the outer parts of the non-Keplerian transition region, and is 
larger than 2 in the inner parts. Fig. 1 shows the run of il(r) 
and n{r) in a typical transition region. The numerical data for 
this example are obtained from the model curve shown in the 
panel (a) of Fig. 2 in Erkut & Alpar (2004). 

For a viscous accretion disk, we write the turbulent viscos- 
ity. 



(41) 



using the a-prescription (Shakura & Sunyaev 1973). The an- 
gular velocity profile in a disk depends on the radial profile of 
the vertically integrated dynamical viscosity. 



t^So = —f(r), 
in 



(42) 



and the appropriate boundary conditions. Here, f(r) is a di- 
mensionless function for the dynamical viscosity. The run 
of / as a function of the radial distance in the typical disk- 
transition region is also shown in Fig. 1 (see Erkut & Alpar 
2004). The radial profile of / is independent of the particular 
prescription for i/. However, once we specify the kinematic 
viscosity as in equation (l4Tl i. we can relate the surface density 
profile (see equationfTTIi. using equations i4T[ and (l42T i. to the 
orbital frequency profile and other parameters of the boundary 
region through 



/3 = 



dlnf dlnn 
■ + ■ 



d\nc. 



dlncs 



(43) 



dlnr d\nr dlnr dlnr 
The first two terms, which we delineate as /3o characterize 
each of the boundary region models developed in Erkut & 
Alpar (2004), with values in the range -15 < /9o < 3. 

Using equations (Q, (1411 1. and (l42l i. we obtain the explicit 
dependence of the ratio of the radial drift velocity and the 
sound speed on the disk-boundary region parameters as 

^1/ _ _VrO _ 3 / a\ f|s , ■ 

n,~cs~2{fjn' 

Note that this ratio, for the given values of a and ilj/ri, takes 
its maximum value for the minimum value of /. As fli,, which 
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determines the growth rates of the modes, depends on the ratio 
rjj/il, the growth rates differ for sufficiently large hydrody- 
namic corrections, i.e., Slj < il, and the fastest growing modes 
can be identified. This is the regime of only slightly super- 
sonic azimuthal flow. In a magnetic boundary layer or transi- 
tion zone where il < JIk, the speed of sound can be written as 
C.5 ~ V,KVHr whereas in the weakly magnetized outer disk, 
Cj ~ r^K^^ (see Erkut & Alpar 2004). We therefore expect to 
find larger values for ils/^ in a non-Keplerian boundary re- 
gion; njn ~ (f7K/r2)(i//r)i/2. The larger the ratio njn, 
the stronger is the effect of hydrodynamic corrections on the 
growth rates of the modes. 

We explore the run of the frequencies and the growth rates 
of free oscillations excited at a characteristic radius r in the 
innermost disk as the hydrodynamic parameter fls/il changes 
at that radius. As a guideline, we have employed the model 
solutions of Erkut & Alpar (2004). These models qualitatively 
represent conditions in a boundary-transition region. The key 
model parameters can be easily translated into the parame- 
ters in equations ( l40b . ( l43T l, and ( l44b . The zeroth-order mag- 
netic force in equation ( l23T l can also be estimated from these 
model solutions with a range < (F™^)o/v,.ori < 1 in the 
inner boundary region. The maximum value of this param- 
eter corresponds to the radius where il is maximum. Near 
the innermost radius of the disk, the azimuthal magnetic force 
vanishes as the relative shear between the magnetosphere and 
the disk becomes negligible, i.e., 17 -fi* ~ 0. 

To elucidate how the hydrodynamic parameters give differ- 
ent growth rates for different modes, we investigate the excita- 
tion of free oscillations in the typical boundary region shown 
in Fig. 1. To see the hydrodynamic effects on the growth 
rates of the modes, we use the full eigenfrequency solutions 
for axisymmetric and non-axisymmetric perturbations given 
in equations ( IA6I )-( IA8I ) and (IA22b - (IA241 i (see Appendix). 
The panel (a) of Fig. 2 shows the real parts of the complex 
mode frequencies at the radius where is maximum for the 
m = and m = 1 cases. For the boundary layer model cho- 
sen, the model parameters at this radius are k/Q. = 2, j3 = -14, 
/ = 0.067, and (F"''°)o/vrof7 = 1. The panel (a) of Fig. 2 is 
obtained for a = 0.1. Although the real parts of the mode fre- 
quencies do not change under a different value of the viscosity 
parameter, say a = 0.01 at this radius, the imaginary parts do 
(see panel b of Fig. 2). The unstable modes grow faster for 
a = 0.1 than they do for a = 0.01, as shown in Fig. 2 (panel 
b). This is because the growth rates are primarily determined 
by ri^/fij oc a// (see equation |44]|. In Fig. 3, we display the 
real (panel a) and imaginary (panel b) parts of the mode fre- 
quencies at the innermost disk radius where = for m = 
and m = 1. The model parameters illustrated in Fig. 3 are 
n/n = 2.5, (3 = -5, / = 1, and (F™^)o/v,ofi = 0. Fig. 3 is 
plotted for a = 0.1. Note that the growth rates, for the same 
value of the viscosity parameter a, are lower than those at the 
radius where is maximum. 

The analysis for the particular transition zone discussed 
above can also be extended to other sample boundary re- 
gions with different model parameters. Our conclusion that 
the fastest growing modes are excited near the radius r = ro 
where ft is maximum remains valid for all sub-Keplerian in- 
ner disks around slow rotators: The common property of all 
boundary regions is that the parameters |/3| and a/f, which 
determine the growth rates of the unstable modes, attain high- 
est values near ro where k/^ ~ 2. Note that the growing 
mode frequencies significantly exceed test-particle frequen- 



cies in the boundary region for ilj/fi < 1 . In the limit of small 
hydrodynamic corrections (rij/fi <C 1), the mode frequencies 
converge to the corresponding test-particle frequencies (see 
panel a of Fig. 2 and panel a of Fig. 3) with negligible growth 
rates as expected (see panel b of Fig. 2 and panel b of Fig. 3). 

The frequency branches that map to the test-particle fre- 
quencies Luf^ = K and Lu[^2 = K±fl in the limit of small hydro- 
dynamic corrections, n jil <C 1, have positive and rather fast 
growth rates at r < ro, where k/Q > 2, when ilj/ri < 1 (see 
panel a of Fig. 2 and panel a of Fig. 3). We identify ujf^ = 

and tjj'* ~ f2 as the quasi-stable or decaying mode frequen- 
cies in the same range of fts/ft as shown in Fig. 2 (panel b) 
and Fig. 3 (panel b). Note that while each of the growing 
branches differs from test-particle frequencies as ils/ft in- 
creases, the difference frequency between consecutive bands 

is always Alu = uj\^'' -wf^ < J7 or Aw = Wj^^'-Wj'' ^ ^ over 
a wide range of fij/ri (see panel a of Fig. 2 and panel a of 
Fig. 3). As ri > ri* in a boundary layer, we find Aui ~ fi, as 
observed for the relatively slow rotators among the neutron- 
star LMXBs that exhibit kHz QPOs. 

4. DISCUSSION AND CONCLUSIONS 

We have studied the global hydrodynamic modes of long 
wavelength in the boundary or transition region of viscous ac- 
cretion disks as a possible source of kHz QPOs in neutron-star 
LMXBs. The stability of the eigenmodes strongly depends on 
the global disk structure imposed by the boundary conditions 
as expected from the mode analysis in the long wavelength 
regime. Our local treatment takes account of the local effects 
of the global disk parameters on the excitation of hydrody- 
namic free oscillations. 

We find that the frequencies and growth rates of the modes 
are mainly determined by global disk parameters such as n/ft, 
the ratio of the radial epicyclic frequency to the orbital fre- 
quency, the radial surface-density profile /3, r^y/fi,, the ratio 
of the radial drift velocity to the sound speed, and ils/ft, the 
ratio of the sound speed to the azimuthal velocity. The local 
values of k/H and /3 directly follow from the global solution 
for the rotational dynamics of a boundary-transition region 
model. The parameters il,y/fts and ils/il depend additionally 
on the particular viscosity prescription and the details of the 
structure of accretion flow in the inner disk. 

The hydrodynamic effects on the frequencies and growth 
rates of the modes are reflected by the parameters Sl^/rij 
and njil. We have found that the growth rates of differ- 
ent modes for a given azimuthal wavenumber m differ sig- 
nificantly when the hydrodynamic corrections are sufficiently 
large, i.e. il^/fts < 1 and ils/ft < 1. In this regime, the 
growing mode frequencies significantly exceed test-particle 
frequencies for a plausible range of fls/ft at any particular ra- 
dius within the boundary region. This is because the effect of 
hydrodynamic corrections on the growth rates of the modes is 
expected to be strong in a magnetic boundary layer or transi- 
tion zone as discussed in § 3.2. In the limit of small hydrody- 
namic corrections, i.e. Ojy/^, <^ 1 and Oj/fi ^ 1, the eigen- 
modes have test-particle frequencies with negligible growth 
rates. Our analysis shows that taking account of hydrody- 
namic effects has an important outcome: the frequencies of 
the growing modes gain higher values above test-particle fre- 
quencies as their growth rates increase for larger hydrody- 
namic corrections. 

Throughout the boundary region, we have identified the 
growing mode frequencies by the eigenvalues that approach 
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the test-particle frequencies = k, ui = n + fl, and = 
K-fi in the limit of small hydrodynamic corrections (see 
panel a of Fig. 2). Modes such as ujf^ = and cj^'' ~ ft are 
not excited for the range of used in Fig. 2 (panel b) and 
Fig. 3 (panel b). The difference frequency between succes- 
sive bands of growing modes is Auj ~ il, as observed for the 
relatively slowly rotating neutron stars in LMXBs that show 
kHz QPOs (§ 3.2). In a boundary layer or transition region, 
the orbital frequency of the disk matter is close to the rotation 
frequency of the star (il ^ 17,). As Auj ~ il, the separation 
between consecutive mode frequencies is always related to 
the spin frequency of the neutron star. This is the result of the 
boundary condition imposed by the rotating magnetosphere 
on the innermost disk. The modes of the test-particle frequen- 
cies Luf' = K and cjj'' = K-fl were originally proposed by Al- 
par & Psaltis (2005) to be associated with the upper and lower 
kHz QPOs observed in neutron-star LMXBs. The frequency 
separation between the higher and lower-frequency kHz QPO 
peaks decreases by a few tens of Hz as both QPO frequen- 
cies increase by hundreds of Hz in a range of 200- 1200 
Hz per source in all observed sources, namely, Sco X-1, 
4U 1608-52, 4U 1728-34 and 4U 1735-44 (van der KHs 
2000). This pair of frequencies satisfies the observed behavior 
that the difference frequency Auj between two kHz QPOs de- 
creases as both frequencies increase. According to our present 
analysis, the pairs of observed kHz QPO bands could be ei- 
ther tjj'* ~ K+ 17 and ujf-' ~ k or uif^ ~ k and Wj'' — K-ft, re- 
spectively, since the identification of the observed QPOs with 
wj'' ~ K + il and ujf^ ~ k also meets the observed condition 
that Auj decreases when both QPO frequencies increase. The 
growth rates of free oscillations produce a wide spectrum of 
sidebands, tjj™' = k± \m\il, some of which fall in the range 
of observed power spectra. The reduction to only two of kHz 
QPO bands as observed in neutron star sources must therefore 
be a product of forced oscillations, resonances, and boundary 
conditions. The role of forced oscillations and resonances in 
the excitation of particular modes will be discussed in future 
work. 

We find that the fastest growing modes with frequency 
branches k and k ± J7 are excited near the disk radius r < ro, 
where the orbital frequency is maximum. In the boundary 
layer, the loss of centrifugal support near ro leads to some ra- 
dial acceleration of the disk matter. The subsequent rise in 
the radial drift velocity of the disk matter is accompanied by a 
sudden drop in the surface density as the vertically integrated 
dynamical viscosity is minimized (see § 3.2). This picture is 
common to all boundary layers or transition zones for which 
magnetic braking is efficient (see Erkut & Alpar 2004 and 
references therein). Thus, independent of the boundary con- 
ditions and the magnetic field configuration in the boundary 
region, there exists a specific disk radius where the hydro- 
dynamical background quantities such as the surface density 
So and the radial drift velocity v,o change dramatically. The 
steepness in the change of the surface density in the radial di- 
rection is reflected through the parameter (3 (see equationfTTTi. 
The sign and magnitude of /3 (the surface density profile in 
the disk) control whether the free oscillation modes will grow 



or decay. We see that the modes of frequency bands k and 
Kztn can be excited with the largest growth rates only for 
the disk radii r < ro, where f3 is minimum and ft is maximum 
(see § 3.2 for the range of /?). The inverse timescale, f2^, 
associated with radial accretion also obtains its highest value 
near the same radii yielding the maximum growth rates in the 
accretion disk. 

The boundary layers with sub-Keplerian rotation rates are 
usually expected to be realized in the innermost regions of ac- 
cretion disks around slow rotators. For a slow rotator, the Ke- 
pler rotation rate at the inner disk radius flKifin) prevails over 
the stellar rotation rate 11*. As we have illustrated in § 3.2, 
the difference frequency between two consecutive bands of 
growing modes in the boundary region is close to the spin fre- 
quency of the neutron star (Aut ~ 11,). This result agrees with 
the kHz QPO observations of the relatively slow rotators with 
spin frequencies below ~ 400 Hz in neutron-star LMXBs. In 
these sources, for example, in 4U 1728-34, the lower kHz 
QPO frequencies increase from 600 to 900 Hz and the upper 
kHz QPO frequencies increase from 950 to 1200 Hz, while 
the difference between the two observed kHz QPO frequency 
bands decreases slightly, of the order of 10 Hz, while remain- 
ing still commensurate with the spin frequency of the neutron 
star (Mendez & van der Klis 1999). A second class of kHz 
QPO sources, such as the sources KS 1731-260, Aql X-1, 
4U 1636-53, which have spin frequencies above ^ 400 Hz 
is usually depicted as fast rotators (see Wijnands et al. 2003). 
In this second class of sources, for example in 4U 1636-53, 
the lower kHz QPO frequencies increase from 900 to 950 Hz 
and the upper kHz QPO frequencies increase from 1150 to 
1 190 Hz, while the difference between the two observed kHz 
QPO frequency bands again decreases by amounts of the or- 
der of 10 Hz (see van der Klis 2000 and references therein). 
The frequency separation between the kHz QPOs is close to 
half the spin frequency of the neutron star if it is a fast ro- 
tator. The emergence of two seemingly different classes of 
neutron-star LMXBs could be related to whether the rotation 
rate of the neutron star is less or greater than that of the inner 
disk matter. It is probable that > HkCHh) for the relatively 
fast rotating neutron star sources for which the difference fre- 
quency between twin QPO peaks is Auj ~ ll,/2. Recent work 
by Mendez & Belloni (2007) actually suggests a more con- 
tinuous distribution between Auj ~ H, and Aut ~ H,/2. The 
rotational dynamics of accretion flow in the innermost regions 
of accretion disks around these fast rotators could be quite dif- 
ferent from those of the transition regions we consider here. 
We plan to address this issue in a subsequent paper. 
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APPENDIX 
GENERAL EXPRESSIONS 

The dispersion relation of the global force-free hydrodynamic modes we consider in § 3 is 

(uj-mfl)[(uj-mflf-fll]-fll-i[(uj-mflfflc + nl] =0, (Al) 
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nl = KH [m^ + /3( 1 + /3)] + [/?( 1 + /3) - ( 1 + 2/3)( 1 - 7.)] , 

.2 



ni = 2m( i+p)nin+ mpn^ 



in 



(A2) 

(A3) 
(A4) 

(A5) 

Note that the solution of equation ( lAlb . in the ideal case of an inviscid fl^, = 0, cold disk ilg ~ 0, where the gas particles barely 
interact with each other, is given by wj™' — i ^"d uj^^^ ~ mil. In general, equation ( lAll i. for the realistic case offl^^/O and 
yo, yields 



and 



= pn^K^ - /3(i + /3)(i - - (i + /3) [/3(i - 7,) - 



,(0) 



.a 



and 



,(0) . 



-III 2- ■ - 3 - 2 

as the eigenfrequency solutions for axisymmetric perturbations (m = 0), where 
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For small hydrodynamic corrections, 772 < »7i, we obtain 

(yTTi7+e2)'^'r!„ (l-£3)^^a 
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as the only non-degenerate eigenfrequency solutions to the dispersion relation for axisymmetric perturbations, where 
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£1 = ^ 7^ + 



2m 



4 \n. 



£2 ^ 



3V3 fnA\V3n^_Vl(qc_y 



and 
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(A16) 
(A17) 
(A18) 

(A19) 
(A20) 

(A21) 
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The general solutions to the dispersion equation for m ^ are 



Abramowicz, M. A., Karas, V., Kluzniak, W., Lee, W. H., & Rebusco, P. 2003, 
PASJ, 55, 467 

Alpar, M. A., Hasinger, G., Shaham, J., & Yancopoulos, S. 1992, A&A, 257, 
627 

Alpar, M. A., & Psaltis, D. 2005, preprint ^astro- ph/0511412^ 
Alpar, M. A., & Yilmaz, A. 1997, NewA, 2, 225 
Erkut, M. H., & Alpar, M. A. 2004, ApJ, 617, 461 
Kato, S. 2001, PASJ, 53, 1 

Mendez, M., & Belloni, T. 2007, MNRAS, 381, 790 
Mendez, M., & van der Klis, M. 1999, ApJ, 517, L51 
Miller, M. C, Lamb, F. K., & Psaltis, D. 1998, ApJ, 508, 791 



(A22) 



(A23) 



(A24) 



(A25) 



(A26) 



(A27) 



(A28) 



(A29) 



(A30) 



(A31) 



(A32) 



(A33) 



(A34) 



Psaltis, D., Belloni, T., & van der Klis, M. 1999, ApJ, 520, 262 
Psaltis, D., & Norman, C. 2000, preprint I astro-ph/000 1391} 
Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337 
Stehle, R., & Spruit, H. C. 1999, MNRAS, 304, 674 
Stella, L., Vietri, M., & Morsink, S. M. 1999, ApJ, 524, L63 
van der Klis, M. 2000, ARA&A, 38, 717 
Wagoner, R. W. 1999, Phys. Rep., 311, 259 

Wijnands, R., van der Klis, M., Homan, J., Chakrabarty, D., Markwardt, C. 
B., & Morgan, E. H. 2003, Nature, 424, 44 



(m) , 



L.J -mQ + (^+A)fli, + i—, 



and 



where 



- - 



A = 



1 /o 



1 



1/3 



2 f nln. 



+ 2 



and 



27 \nh I ' n \ nl 



_^ 4 (^aW 1 (^^^'^ 



27 [vii, 



In the limit of small hydrodynamic corrections, \C,i \ ^ we find 
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as the only non-degenerate eigenfrequency solutions to the dispersion relation for non-axisymmetric perturbations. 
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Fig. 2. — The real and imaginary parts of the solutions to the dispersion relation for axisymmetric (m = 0) and non-axis^inmetric (m = 1) perturbations excited 
at the radius where the orbital frequency Q is maximum. The model parameters are k/H = 2, /3 = —14, / = 0.067, and (f^ ''^)o/vrO^ = 1- The curves are labelled 
with the values of the frequencies in the test-particle limit. The k and ft ± 11 modes have the same positive growth rates, growing more rapidly for larger a. The 
zero-frequency and SI modes, indicated by the overlapping dotted curves, do not grow for any a. 
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Fig. 3. — The real and imaginary parts of the solutions to the dispersion relation for axisymmetric (m = 0) and non-axisymmetric (m = 1) perturbations excited 
at the innermost disk radius where f2 = Q». The model parameters are k/H = 2.5, /3 = -5, / = 1, and (F^''*)o/vH)n = 0. The frequencies and growth rates are 
obtained for a = 0.1. The growing modes are associated with the K-t-H, k, and k-H branches, with k + H having the highest growth rate. The zero-frequency 
and f2 modes grow for lower values of fij/n, but decay in the regime of fJj/Q < 1, where the k and k ± f2 modes have the fastest growth rates. 



